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ABSTRACT 

We find the forms of the transfer equations for polarized cyclotron radiation in the 
atmospheres of compact stars, which are simple enough to allow practical implementa¬ 
tion and still preserve all important physical effects. We take into account a frequency 
redistribution of radiation within the cyclotron line as well as the relativistic and 
quantum-electrodynamic effects. Our analysis is valid for the magnetic fields up to 
10 13 G and for temperatures well below 500 keV. We present and compare two forms 
of the radiation transfer equations. The first form, for the intensities of ordinary and 
extraordinary modes, is applicable for the compact stars with a moderate magnetic 
field strength up to 10 11 G for typical neutron star and up to 10 9 G for magnetic white 
dwarfs. The second form, for the Stokes parameters, is more complex, but applicable 
even if a linear mode coupling takes place somewhere in the scattering-dominated at¬ 
mosphere. Analysing dispersion properties of a magnetized plasma in the latter case, 
we describe a range of parameters where the linear mode coupling is possible and 
essential. 

Key words: line: formation - radiative transfer - scattering - stars: neutron - white 
dwarfs. 


1 INTRODUCTION 

An emission from the atmospheres of compact stars (neutron stars and white dwarfs) can provide a valuable information on 
the physical properties of plasma surrounding these objects. The most universal feauture of that emission is the cyclotron 
line, whose formation was studied by many authors (e.g., Pavlov et al. (1980); Wang et al. (1989); Bezchastnov & Pavlov 
(1991); Araya-Gochez & Harding (2000)). A significant progress in the X-ray observations, facilitated by the space-based 
observatories like XMM-Newton and Chandra, revives an interest in the theoretical and observational studies of a line formation 
in the atmospheres of such stars (Sanwal et al. 2002; Gotthelf & Halpern 2009; Pavlov et al. 2004; Mushtukov et al. 2015b). 
A discovery of the spectral features in the thermal spectra of the X-ray dim isolated neutron stars and in the central compact 
objects (CCOs) is especially interesting. One of the possibilities is that they are the cyclotron lines forming in the outer layers 
of the atmospheres (see, e.g., Potekhin (2014); Marsden et al. (2001); Triimper et al. (2013); Haberl et al. (2004)). 

The physical parameters in the upper (near-photosphere) layers of the compact stars’ atmospheres are such that scattering 
of photons strongly dominates over their absorption (Wang et al. 1989; Zheleznyakov 1983), that is an electron collision rate is 
much smaller than a spontaneous radiative rate. In this case, a population of the Landau levels is governed by a local number 
density of the resonant photons. Usually, an occupation number of resonant photons in the upper atmosphere is much smaller 
than a thermodynamic equilibrium value, so that the excited Landau levels are underpopulated to such an extent that almost 
all electrons reside at the ground level and an effective transverse temperature is nearly zero (Nagel & Ventura 1983). The 
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longitudinal momentum distribution of electrons is supposed to be Maxwellian (in line with most of authors), though in the 
upper layers of the atmosphere, where the plasma is almost collisionless, this assumption has to be taken with care. 

A strong magnetic field of the neutron stars makes the quantum-electrodynamic (QED) effects also important (Pavlov & Shibanov 
1979). A number of authors investigated in detail an influence of the QED effects on the dispersive characteristics of electromag¬ 
netic waves in a magnetized plasma (see, e.g., Pavlov & Shibanov (1979); Adler (1971); Kirk & Cramer (1985)). It was shown 
that a vacuum polarization can affect the radiation transfer by modifying a polarization of the normal waves (modes) and induc¬ 
ing the additional features in their opacity coefficients (Ho & Lai 2001; Zheleznyakov et al. 1983; Zhelezniakov & Litvinchuk 
1983). Also, it was indicated in Lai & Ho (2003), that, due to a linear mode coupling, a description of a radiative transport 
within the framework of geometrical optics (which is widely used in lots of calculations of a radiation transfer in the atmo¬ 
spheres of compact stars) is, in general, inadequate near the points, where the vacuum contribution to the refraction indices 
becomes comparable with the plasma one, i.e., in the vicinity of the so-called ’vacuum resonance’ points. Below we show 
that, even in a relatively weak (for the neutron stars) magnetic field of the order of B ~ 10 11 G, the vacuum polarization can 
significantly modify a cyclotron line profile in the emergent spectra. 

A major part of a prior work on the modeling of cyclotron lines in the spectra of highly magnetized compact stars has 
been done either within a cold plasma approximation (Shibanov et al. 1992; Suleimanov et al. 2009, 2010) or within a limit of 
a nonrelativistic cyclotron resonance condition (Bespalov & Zheleznyakov 1990; Zheleznyakov et al. 1999; Suleimanov et al. 

2012). The latter approach greatly simplifies a solution since the radiation transfer equations can be reduced to the quasi- 
monochromatic ones by a proper choice of variables. As it was firstly mentioned in Wasserman & Salpeter (1980), at each 
resonant scattering there is a conserved quantity - a velocity of an electron, which is in resonance with a given photon, 


P = 


UJ — Wb 


(1) 


oj cos 6 ’ 

where w is the frequency of a photon, 9 the angle between the photon wave vector and the magnetic field, ujb — eB/m^c 
the cyclotron frequency, and /3 the velocity of electron in the units of the speed of light c. This allows one to simplify a 
radiation transfer problem by using a velocity (3 instead of a frequency u> as an independent variable. Such an approach proved 
to be very helpful in the limit of low optical depths and was widely used for studying the cyclotron line formation in the 
neutron star magnetospheres (Lyutikov & Gavriil 2006; Rea et al. 2008). In what follows, we refer to that approximation as a 
quasi-coherent approximation, because the changes in frequencies due to a Doppler effect do not accumulate and all resonant 
photons remain locked in the line core. 

Recently we showed (Garasyov et al. 2008, 2011b) that, for the extremely magnetized white dwarfs and for the neutron 
stars with magnetic fields stronger than 10 10 G, the relativistic corrections to the Doppler resonant conditions become impor¬ 
tant and brake the quasi-coherent approximation, leading to an efficient escape of the cyclotron radiation from a relatively 
large optical depth. In their most general form, the transfer equations for a magnetized plasma + vacuum medium (see, e.g., 
Mushtukov et al. (2012)) are too complicated even for numerical simulations. 

The goal of this paper is to find the simplest possible forms of transfer equations, which would be simple enough for 
practical implementation and, at the same time, still include all the relevant physical effects and provide accurate results. 
The minimal set of physical effects to be incorporated into the transfer equations depends on the physical conditions in the 
atmosphere of compact stars. We explore the parameter space for a fully ionized magnetized plasma in a mildly-relativistic 
approximation (with temperature in energy units T <C m e c 2 ) with moderate (up to 10 13 G) magnetic fields, and identify which 
form of the transfer equations is suitable for which parameters. Our formulation of the radiation transfer problem allows for 
the effects of the vacuum polarization as well as for a non-trivial redistribution of photons over frequencies and their escape 
from the cyclotron line due to the relativistic Doppler effect in the scattering-dominated atmospheres. 

We concentrate on describing correctly the scattering of photons near the fundamental harmonic of cyclotron resonance. 
This harmonic is the most prominent feature in the spectrum, and radiation pressure at the fundamental frequency is a 
major factor (Mitrofanov & Pavlov 1982; Mushtukov et al. 2015a; Wielgus et al. 2015), which determines the structure of the 
atmospheres and possible formation of radiation-driven winds. For this reason, we do not consider higher harmonics, but 
adding them is straightforward. 

The paper is organized as follows. In the next section, we derive all the coefficients in a complete set of the radiation 
transfer equations for the incoherent mode intensities. The coefficients are expressed through the dielectric tensor for the 
’warm mildly-relativistic plasma + magnetized vacuum’ media, which is also calculated. That set of equations can be used 
only in the case of the geometrical optics, where the radiation transfer in an inhomogeneous plasma is adiabatic, i.e., preserves 
a wave’s mode composition in the absence of scattering. In the third section, we derive a more general but more complex set 
of the transfer equations for the Stokes parameters. They are to be used when the adiabatic approximation is not valid. In the 
last section, we discuss the range of plasma parameters, where the radiation transfer can be non-adiabatic, i.e., can undergo 
the linear mode coupling. We show that the vacuum resonance zone for the cyclotron radiation is located at the higher optical 
depths as compared to the continuum and may lie far below the photosphere. We also indicate the three domains in the 
parameter space, related to the three approximations useful in the radiation transfer problem. In the discussion section, we 
outline the most important results of the work and their possible astrophysical implications. 
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2 TRANSFER EQUATIONS FOR THE X-, O-MODE INTENSITIES 


The simplest, yet widely employed way to describe the transfer of a polarized radiation in the magnetized plasma is to 
use the transfer equations for the intensities of two electromagnetic modes - ordinary(O) and extraordinary (X) ones. In a 
tenuous medium, where a difference of the dielectric tensor from its vacuum value is small, any electromagnetic wave is almost 
transverse and can be easily represented as a superposition of two modes. The magnetoactive plasma is a tenuous medium at 
frequencies well above the plasma frequency oil = \/4:ne 2 N e /m e (and outside of strong resonances), where e is the elementary 
charge, N e the number density of electrons, and m e the mass of an electron. In this case the scattering phenomena became 
transparent if described in the form of two incoherent modes. The extraordinary mode has a resonance in opacity in the 
vicinity of cyclotron harmonics. The ordinary mode also has a resonance, but it is much weaker and completely disappears 
in the limit of a cold plasma. In the geometrical optics (adiabatic) approximation, these modes are coupled with each other 
only through scattering. 

The radiation transfer equations in this case are (Nagel 1981) 


dJi , (s 
d7 = - (<7 ‘ 


:) + af bs) )Ji + 


E 

3=1 


da. 


(sc) 


u 1 dfldtu 


(w',n' 


w,n)Jj(u/,n') dO'du/ + ^R,,T<rf bs) . 


( 2 ) 


Here Ji(Q,io,s) is the intensity of radiation in the i-th mode with a unit polarization vector ei related to a unit spectral 
interval and a unit solid angle (a usual convention is i — 1 for the extraordinary mode and i = 2 for the ordinary mode), 
dfi the elementary solid angle, s the coordinate along the ray, dcrj^ /dfldcn the differential scattering cross section, which 
describes scattering from mode j to mode i, cr| abs ^ the coefficient of the free-free absorption for the i-th mode, af^ = 
f 52j =i 2 do»j/(dfIdw)df2dw the total scattering cross-section, and 

d = _ 1 / o\ 

U1 ' T 47r 3 c 2 exp (hw/T) — 1 

the thermal equilibrium (Planck) intensity. The primed quantities are for the incident wave, those without prime are for the 
scattered wave. Here and below T is a temperature measured in energy units. 

A photon redistribution over frequencies in the cyclotron line is, in general, an important effect (Garasyov et al. 2011b), 
which can lead to a significant (by several orders of magnitude) increase in the escape probability for the resonant radiation 
from large optical depths. The frequency redistribution in the cyclotron line is due to the relativistic Doppler effect. A motion 
of electrons in the strong magnetic field is essentially unidirectional, and, hence, the electrons resonantly interact only with 
the photons of almost certain frequency. The resonance condition for scattering by an electron, which undergoes a transition 
between the ground level and the £-th Landau level, takes the form (Kirk & Cramer 1985) 

huj 


m e c‘ 


+ s/l + P 2 - y / l + (P + fifc||) 2 + 26£ = 0, 


(4) 


where P is the momentum of an electron along the magnetic field (in the units of m e c), k\\ the component of the photon’s 
wave vector along the magnetic field direction (in the units of m e c/h), b = B/B CT the magnetic field strength (in the units 
of Schwinger field B CI = mfc 3 /(he)). Eq.(4) has two roots, whereas the non-relativistic resonance condition (1) has only one. 
The use of only one root, instead of two, does not allow for the actual redistribution of photons over frequencies, i.e., for a 
change in a resonant group of electrons, and prevents from their escape from the cyclotron line core. This leads to inaccurate 
results even if the plasma temperature is far below the relativistic values. 

The simplest option to preserve the two roots in the resonance condition is to expand Eq. (4) into a power-law series, 
keeping terms up to the second order in /3: 


w(l — 0 cos 9 — -f—b cos 2 9) = £ujb( 1 
2 lob 



(5) 


The resonance condition in such a form describes all essential effects near the cyclotron resonance in a mildly-relativistic 
plasma. In what follows in this section we present the scattering properties and the wave characteristics in a warm plasma. 
Let us assume that the electron distribution function over the longitudinal velocities is Maxwellian, 

7)2 


m = 




■ exp < - 


20 2 


( 6 ) 


where 0t = \/T/(m e c 2 ) is a typical thermal velocity of electrons (in the units of the speed of light). Also we assume that 
all the electrons are at the ground Landau level and consider only excitations to the first Landau level, that is l = 1 in 
Eq. (5). Usually it is safe to neglect the population of excited Landau levels in the upper layers of compact stars’ atmospheres 
(Nagel & Ventura 1983). The reason is that the rate of spontaneous emission is many orders of magnitude higher than the 
rate of collisional excitation. 

To derive the transfer equations, which take into account the frequency redistribution of photons, let us start from the 
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general expressions for the differential scattering cross-section (Ventura 1979) 

d(j( SC ) / r 2 

dflcL; = r o~^ J |( e 'i|n|e;)| 6(u - to' + Aio), 


dfldto 

where ro = e 2 /(m e c 2 ) and 
Ato = —to/3 cos 9 + to' f3 cos 6' — b 


to a to a’ 

— cos 9 -cos 9 

LOb OJb 


It is also convenient to express the total scattering cross-section in the same notation, 


o-f c) = Y r o / |fie 

0U£ 

= —47rro—Im(e*|n|ei). 


8n 


I 


and through the optical theorem 

>c) = _ 4nn C 
tO 

In the dipole approximation, the matrix elements of the scattering amplitude for an individual electron is equal to 


to 


iixx — riuy — 


l 


l 


2 Vw + cub to* — tot 


CU + C Jb 


II XZ — Ilza; — in 1J Z — in??/ — 


to f to/3 sin 9 


U zz = 


to 


to + iq 
to* = w(l — (3 cos 9 — 


hto 

2 m e c 2 


2cob \to* — tot 


cos 2 9) + iq, 


tOg — Wb(1 — f3 2 /2 — -). 


(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 


These equations are derived as a mildly-relativistic limit from the general form which can be found in Meszaros (1992) 
(equation 5.3.5). We use the reference frame with z-axis along the magnetic field and x-axis in the plane, containing the 
magnetic field and the wave vector (the same reference frame will be used below unless specified otherwise). In the above 
expressions, 7 = q r + 7 fF is the total resonance width for an individual particle, which takes into account both the radiative, 
7 r, and collisional, 7 ff, broadening 1 . Based on the fact that these partial line widths are related to the total scattering cross 
section and to the free-free absorption, respectively, Pavlov & Kaminker (1975) deduced that the ratio between the scattering 
cross section and the absorption cross section equals to the ratio of the corresponding line widths. Thus, the cross section for 
the free-free absorption is 


(aba) 7ff (sc) 
a i = - a i ■ 

7r 


(17) 


In the non-relativistic limit, j r = (2/3)e 2 u; 2 /(m e c 3 ), where «f = e 2 /(fic) is the fine structure constant. In the quantum 
approach, which is necessary to use in the stronger magnetic fields, the radiative damping rates was calculated in Gonthier et al. 
(2014). A detailed evaluation of 7 ff can be found in (Potekhin & Chabrier 2003); in the most interesting frequency range near 
the cyclotron frequency, 


2n at to\ 

7 * ~ V 27 — 


'- exp i-£) 


where u = T/(hto), and g a is the Gaunt factor. In the limit of htos <C kT one has 


9 -^r xp 


(— 

\2hto 


K 0 


( — ) 

\2huj J ’ 


(18) 


(19) 


where Kq is the modified Bessel function, and the expressions for arbitrary magnetic fields could be found in (Potekhin & Chabrier 
2003). 

The polarizations of the normal modes can be derived from the dielectric tensor of plasma + vacuum media e, which is 
related to the matrix II (Meszaros 1992): 

e = i + x (vac) +x (pl) ; (20) 


1 A more refined, but still approximate, methodology calls for the use of two different collisional frequencies for the longitudinal and 
transverse motions (Potekhin Sz Chabrier 2003). Since this difference can hardly lead to an observable effect in the scattering-dominated 
atmospheres, we adhere to a simpler approach. 
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x 


(pi) _ 


_^L f 


n/(/3)d/?; 


( 21 ) 


x( vac ) = —2a(k ® k) + 4a(kxB) ® (kxB) + 7a(B ® B), (22) 

where a = af6 2 /(457r). The expression (22) for permittivity tensor of vacuum is valid for B -C B cr . Integrating Eq. (21) with 
the Maxwellian distribution function (6), we obtain the resulting expression for the permittivity tensor of the magnetized 
plasma 2 

= K b - + b +). 


X ( Tf = ~X=b--b + , 

Y (pi) = y (pi) = _i Y (p!) = i Y (p!) = c 
A.xz A.zx Ai/2 ± A.zy ^— ? 


(pl) = 
Azz 

where 
b+ = 


u>(u} + iy) 


iw 2 


+ d— , 


b _ = 


2lo(uj + vy + ujb) ’ 

v^l ro(Ci) - 
2 w 2 /3| 


£2 - Ci 


f A^u; 2 sin6» Ciro(Ci) - 6^(6) 


2 ujujb^t 


6 - Ci 


d_ = 


sin 2 0 ^ Ci(l + iy^Ci^te)) - &(! + h/tr^rofe)) 

C2 - Ci 


+ i\ 


7i b zu( Ci) — tn(C 2 ) 


16 / 3 : 


C2 - Ci 


(23) 

(24) 

(25) 

(26) 
(27) 


Ci,2 = 


V2/3t 

The function 

+ OO 

™( z ) = ^ J C 


cos 9 ± \l cos 2 6 — 2 I 1 — —- + 


hut 


oj 2 m e c 2 


sin 2 9 + i— 




dy 


(28) 


(29) 


in the region Im(Z) ^ 0 equals to the Kramp function (the plasma dispersion function), whereas for Im (Z) < 0 one can use 
the relation vj{Z') = —vj*(Z*). This function can be easily evaluated numerically with high precision (Schreier & Kohlert 
2008; Letchworth & Benner 2007; Abrarov et al. 2010). Because of the two groups of electrons, which are in resonance with 
a given photon, the resonant terms in the dielectric permittivity tensor in the mildly-relativistic approximation contain two 
plasma-dispersion functions, instead of one in the commonly-used tensor for the warm non-relativistic plasma. 

Introducing a polarization parameter, 


_ —txx cos 2 9 + e yy + sin(20) — e zz sin 2 9 

2i(e xy cos 9 + t yz sin 9) 
and the polarization coefficients of modes, 


Ax,o = q T V<1 2 + 1, 

one can find the polarization of electric field in the X- and O-modes (Meszaros 1992): 

( x,o ) _ Ax,o cos 6 

e * ~ Vl + l^x,o| 2 ’ 


(30) 


(31) 


P (x,o) _ 


(32) 


V Vl + l^x,o| 2 ’ 

(x,o) _ Ax,o sin6> 

\J\ + |Ax,o| 2 ' 

Here the minus sign in the value A corresponds to the extraordinary (X) mode and plus - to the ordinary (O) one. These 
expressions for the polarization vectors are given in the coordinate system where the x-axis is associated with the transverse 
component of the wave vector. This is convenient in the case when the magnetic field is perpendicular to the star’s surface. 


2 MATLAB code, which evaluates the components of the permittivity tensor and refraction indices for a mildly-relativistic plasma, can 
be found at http://www.mathworks.com/matlabcentral/fileexchange/52024 
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Figure 1. Scattering cross sections separately for X-mode (left panel) and O-mode (right panel) for B = 2.2 • 10 12 G, N e = 10 20 cm 3 
and T = 20keV. 


In the general case of the inclined magnetic field, it is more convenient to use another reference frame (let us call it x*y*z*), 
where z*-axis is parallel to the magnetic field and x*-axis lies in the plane that contains both the magnetic field and the 
normal vector to the surface. In that case a wave vector is rotated by some angle from the x*z*-plane. By expressing the 
wave vector k through the polar angles 0 and 0, we obtain the following expressions for the mode polarization components: 


e 

e 

e 


(x,o) 

x* 


( x,o) 

V* 


(X,0) 

z* 


— e^'°^ cos rj> + e y X ’° sin rj>, 

= sin 4> + e| \f'°^ cos <f >, 

- JX.O) 
c 2 


(33) 


We would like to emphasize the importance of xz-, yz-components of the dielectric tensor which are sometimes ignored. 
These components are responsible for the resonance in the ordinary mode. Their omission leads to incorrect opacity for the 
ordinary mode in the cyclotron line core for all propagation angles. 

To compare our results with the results obtained for fully-relativistic case (see, e.g., Araya & Harding (1999)), we calculate 
the scattering cross-sections in the frequency range close to the frequency of the first cyclotron harmonic. It is calculated using 
Eq. (10) for T = 5 — 20keV and plotted in Fig. 1-2. The shape of the resonance for close to transverse propagation angles is 
distorted and acquires a narrow feature with the width of the order of yy due to the presence of a cutoff frequency in the 
resonance condition (5). Near this point the amount of resonant particles per unit frequency interval tends to infinity if we 
neglect the line broadening. The asymmetry in the scattering cross-section is especially pronounced for small values of cos#, 
showing sharper cutoff above the resonance than below it. Note that even at temperatures T <C m e c 2 , the relativistic effects 
change significantly the scattering cross-section of X-, O-modes propagating nearly transverse to the magnetic field lines at 
angles 1 6 — 7t/2| < fir ■ 


3 TRANSFER EQUATION FOR THE STOKES PARAMETERS 

The two equations for the incoherent intensities of the X-, O-modes provide an adequate description of the radiation trans¬ 
fer only if the polarizations of propagating waves evolve adiabatically, i.e., in the absence of the linear mode coupling 
(Zheleznyakov et al. 1983). Generally speaking, the geometrical optics (adiabatic) conditions can break at some point on 
the way out of the atmosphere, and the X-, O-modes become coherent due to the linear mode coupling. In such a case, the 
transfer equations for modes should be replaced by a complete set of the four transfer equations for the Stokes parameters. 

In what follows, we derive the coefficients in this set of equations in the case of the mildly-relativistic highly magnetized 
plasma. Although it is possible to obtain the equations using the matrix n, here we use another technique. For the sake of 
simplicity, we closely follow the approach of Lai & Ho (2003), where the transfer equations were derived for the cold plasma 
and without taking into account the frequency redistribution of photons. 

In the adopted approach, the propagating electromagnetic waves are described by a vector J, which is defined in terms 
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Figure 2. Close-up view of the average total scattering cross section near the cyclotron frequency for b = 0.04 and T = 5keV. The 
narrowest curve corresponds to the almost transverse propagation angle cos# = 0.005. Note the sharp cutoff of the resonance curves due 
to relativistic effects in the case of small values of cos 6. 

of the four Stokes parameters (Meszaros 1992): 

( I + Q \ 


J = 


V 


I-Q 
2 U 
2V 

An evolution of this vector is governed by the transfer equations 

^ = —M3 + S em + S sc , 
as 


(34) 


(35) 


where s is the coordinate along a ray. The source functions S em and S sc are the emissivity and photon yields due to scattering 
from other directions and rays, respectively. The transfer matrix M describes the absorption and the escape of radiation due 
to the scattering as well as the evolution of polarization. The eigenvectors of the transfer matrix correspond to the X- and 
O-modes. 

The components of the transfer matrix M are related to the components of the dielectric tensor of a media: 

21m cn 0 Im?i 2 —Re?i 2 

1 + 7ff/7r o 21m ?22 —Im?i 2 -Re?i 2 

- 2 Im?i 2 2 Im?i 2 Im(?n+<; 22 ) Rc(?ii-c 2 2 ) 

- 2 Re?i 2 -2Re?i2 Re(? 2 2-?n) Im(cn+? 2 2) 


M = 


(36) 


where 

Cn = tax cos 2 9 + e zz sin 2 9 — t xz sin 29, 


( 37 ) 


?22 
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Note the terms t xz and e yz , which are absent in the cold plasma limit. 

Now, let us consider an incident electromagnetic wave with the wave vector k', frequency u>' and electric field Ei n (r, t) = 
Eoe lk r_1 “ t , where Eo = ( Eo x cos 6' cos <f>' — Eo y sin 0 ')xq + (Eq x cos O' sin <\> + Eo y cos </>')yo + Eo x sin 6 'zj. Here Eq x and Eo y 
are the components of the electric field in the reference frame x’y’z’, where the z’-axis is directed along the wave vector k' 
and x’-axis lies in the plane, containing both k' and B. A scattered wave, propagating at angles 9 and <j>, is characterized by 
the wave vector k and frequency lo. In the dipole approximation, the electric field of the wave scattered by an electron is 


E sc = _ [(e x cos 9 cos <t> — E y sin 0)xjj + (E x cos 9 sin <f> + E y cos <j>) yj) + E x sin 9zq] e lkr ,UJt , (38) 

where d is the dipole momentum of an oscillating electron, r the distance between an observation point and an electron. 
Then, for an electron with an undisturbed velocity /3 along the magnetic field, we can construct the matrix, which relates the 
electric field of the scattered wave to that of the incident wave: 

(S) 


where 


ao = A(A0) cos0(cosQ' — /3) + (1 — i—) sin 9 sin 0\ 

uj 

(40) 

bo — T(A0) cos0(l — /3cos9'), 

(41) 

co = T(A0)(/3 — cos O'), 

(42) 

do = A(A0)(1 — p cos 9'), 

(43) 

... o/(l — 3 cos 9') , „ . , . * . . 

A(A 4>) — ---^ 2 —( w cos A</> + io; s sin A<f>), 

UJ* — UJ b 

(44) 

„, . ,. uj(1 — B cos 9') , * . . , .* . ,. 

T(A^) = ——5 - 12 J (u sm A0 — ia; s cos A()>), 

UJ* — UJ R 

(45) 


and A rf> = <f> — <j>'. 

The scattering source function S sc can be expressed in the following way 


a o b o 
Co do 


Eo x 

E 0y 


(39) 


Ssc =rlN e J 

d Cl' J du j dBf{B)R{V -> 

k)J(k , )5(cc — uo' + Aui), 

where the scattering matrix is 





l«o| 2 

M 2 

Re (aobo) 

Im(a3&o) 

R( k' ->• k) = 

M 2 

2Re (a(Sco) 

Idol 2 

2Re (bod 0 ) 

Re (c^do) 

Re (aj^do + 5qCo) 

Im (cgdo) 

Im (aod 0 — feoCo) 


—21m (aoCo) 

—21m (bod 0 ) 

—Im (aodo + 6qCo) 

—Re (aod 0 — &o c o) 


(46) 


(47) 


The emission source function is determined by the requirement that in the thermodynamic equilibrium, when the radiation 
is unpolarized and its total intensity is equal to the Planck intensity, i.e., J eq = (B^,t/ 2, 0, 0), the absorption term 

should cancel out with the term S em - Thus, one finds that 


S 


em 


7fffco-B u,t 
2 7r 


^ Im(oi) 
Im(?22) 
0 

\ — 2Re?i2 


\ 

/ 


(48) 


The equations, derived in this section, are applicable to the transfer of polarized radiation even under the conditions, where 
the geometrical optics of the X- and O-modes is not valid and they are coupled due to propagation in an inhomogeneous plasma. 
The scattering matrix, which we find here, is the simplest possible, which still takes into account all important physical effects. 
It is known that direct solution to the Stokes equations sometimes requires uncomfortably high spatial resolution to produce 
accurate results (Lai & Ho 2003). The main reason for that is the oscillations of V and U Stokes parameters in the vicinity 
of so-called vacuum resonance point in the presence of linear mode coupling, i.e., where the geometrical optics approximation 
fails. A method to overcome this difficulty was suggested in van Adelsberg & Lai (2006), where transfer equations for normal 
modes are separately solved in two regions before and after the coupling layer and joined together by introducing the ’’jump 
conditions”. This method implies no scattering in the coupling layer. However, as we will show in the next section, the 
scattering optical depth of the coupling region (under typical conditions for neutron stars) is greater than or of the order of 
unity for the photons at the cyclotron resonance, so that one has to solve the transfer equations through the coupling layer. 
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4 LINEAR COUPLING OF THE ORDINARY AND EXTRAORDINARY MODES 


The electromagnetic modes cannot be treated independently if their polarization changes too quickly along a ray. Usually this 
happens near the point, where the plasma and polarized vacuum contributions to refraction indices are of the same order. 
This region of so-called vacuum resonance is defined by the condition \q\ 2 < 1. There, the rapid change of polarization vectors 
(32) in inhomogenous plasma can cause non-adiabatic transformation of modes, that is, linear mode coupling (Lai & Ho 2003; 
Zheleznyakov et al. 1983; Zhelezniakov & Litvinchuk 1983). 

In a mildly-relativistic plasma the polarization parameter is 

_ i(b- + b + ) sin 2 9 — 3a sin 2 9 + c_ sin(2f?) — (d_ — oi\/(ijo{uj + if))) sin 2 9 , . 

^ 2i [(&_ — 6+) cos 9 + ic- sin 9] 

and can be simplified in the vicinity of the cyclotron frequency as 
^ sin 2 9 ib_ — 3a 
^ 2 cos 9 i b- 

since the terms c_,d_, &+, and a> 2 /u; 2 are much less than the term Here b- is due to plasma contribution and 3a is due 
to vacuum contribution. Thus, the necessary condition for the mode coupling in the vicinity of the cyclotron resonance is 
16_ | ~ 3a. For estimation purposes, we consider the case, when | cos0| ~ 1. For such propagation angles 

(iO-UJB+Vy\ , . 

8 UJ 2 /3t ™ { y/2UJpT J ' ( } 

Under certain conditions the coupling may be important in the optically thick region. Then the radiation transfer is 
affected by the mode coupling and the use of Stokes parameters is necessary. Let us introduce the scattering optical depth for 
an X-mode photon with a frequency u in a magnetized plasma, 


(50) 


t(ui) = / a x N e (z')dz', 


(52) 


where z is the vertical coordinate in the atmosphere and the axis is directed from infinity to the surface. For order of magnitude 
estimates let us use the simple nonrelativistic expression for scattering cross-section in a pure plasma (without vacuum effects; 
see, for example, Meszaros (1992); Zheleznyakov (1996)) 


ax 


7r 47re 2 (l + cos 2 9) 


Retu 


a ; — ujb + ry 


2 mccj/3 T cos 9 V V2ujf3 T cos 9 

Then, substituting | cos 9 | by unity, we obtain 


r(u}) 


-H^L 

2 c/?tw 


Retxj 


ui — lob + by 


(53) 


(54) 


%/2cj/3t 

where H is the characteristic atmosphere height for an exponential atmosphere with density N ~ exp (z/H). 

As a result, equating the plasma and vacuum contributions, expressing lol via r from Eq. (54), and substituting it into 
Eq. (51), we find the optical depth at which the polarization of the X-mode undergoes transition to the vacuum polarization 
and, hence, the geometrical optics approximation can be broken: 


T, 


12ujHoi{ 
457rc 


B_ 

B CI 




(55) 


The largest value of the function r* (w) is at the cyclotron resonance. If this value exceeds unity, then region of the linear 
mode coupling lies below the photosphere. For example, in an isothermal hydrogen atmosphere with the Boltzmann density 
profile one has 


t*(wb) « 0.2 


T 


( B 


(56) 


v 50eV / V 10 11 G / V 9 

near the resonance (here g is the gravitational acceleration). Thus, for the cyclotron radiation the coupling layer (with possibly 
non-adiabatic polarization change) can be located in the optically thick regions of atmosphere even in the magnetic fields as 
low as 10 11 G, leading to observable effects. At the same time, the linear mode coupling is likely to have a negligible influence 
on the transfer of continuum photons. 

Another important phenomenon which affects the radiation transfer in a resonance line is the frequency redistribu¬ 
tion during scattering. Regarding cyclotron lines, it was discussed in a number of works (Wasserman & Salpeter 1980; 
Zheleznyakov et al. 1999; Garasyov et al. 2011a). Here we briefly estimate the atmospheric conditions under which it will 
play crucial role. Consider a photon, whose probability of escaping from the line due to single excursion into its wing is given 
by Eq. (7) from Garasyov et al. (2011b): 

n Pt Pt 

i esc I - ~ , 

t\J 8 In r 4 r 


(57) 
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Figure 3. The zones of influence for different radiation transfer effects in the parameter space Zone I: the photons leave the 

atmosphere without significant frequency redistribution; in this situation it is sufficient to solve the transfer equations in the quasi- 
coherent approximation. Zone II: the frequency redistribution plays an important role and determines the cyclotron line profile. Zone 
III: additionally, the vacuum polarization is strong enough to change the dispersion properties of the X-, O-modes in the regions, where 
scattering is still important (i.e., below the photosphere, t*(ujb) > !)• The panel a: the neutron stars with g = 10 14 cms -2 ; the panel 
b: the white dwarfs with g = 10 8 cms -2 . 


For the frequency redistribution to become important, this probability should be greater than probability of escaping due to 
spatial diffusion alone, 

Peso > T (58) 

what gives an inequality 


r 



(59) 


This optical depth is to be compared with the so-called thermalization depth (Mihalas 1978), where most of the escaping 
photons are generated, 



implying 


(60) 


b > 5 • 10- 


10 14 cms -2 m e c 2 

a T~' 


(61) 


Figure 3 shows the three zones in the magnetic field - temperature plane, which call for the different approaches to the 
radiation transfer problem. Going from lower to higher zone requires taking into account a larger number of effects and solving 
more complex radiation transfer equations. In zone I, the frequency redistribution of radiation within the cyclotron line core 
is weak, so that the major fraction of emerging photons leaves the atmosphere due to usual diffusion. The radiation transfer in 
this situation can be considered in the quasi-coherent approximation. In zone II, the frequency redistribution of radiation in 
the line’s core is strong, and photons primarily leave the atmosphere due to excursions to the line’s wings rather than due to 
the spacial diffusion. This zone is bound by the condition (61), which ensures that most of the resonant photons in emergent 
spectra have at least one excursion to the line’s wings during their passage through the atmosphere. In this zone one should 
use more complicated transfer equations for incoherent X-, O-modes with partial frequency redistribution effects taken into 
account. Additionally, in the zone III, the vacuum polarization effects become important in the inner regions of atmospheres 
(i.e., where t*(ojb) > 1) and may lead to the linear mode coupling. In this case, the transition between the plasma-determined 
polarization and the vacuum-determined polarization occurs in the optically thick regions of the atmosphere, forcing one to 
use the complete set of transfer equations which can deal with non-adiabatic propagation of modes, i.e., the equations for the 
Stokes parameters. 
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5 CONCLUSIONS 

We analyse, for a wide range of temperatures and magnetic field strengths, the impact of the mildly-relativistic and QED 
effects on the radiative transfer of the gyro-resonant radiation in the atmospheres of compact stars. We find that the parameter 
space can be split into three zones with different requirements to the amount of details in the radiation transfer modeling. In 
the first zone the simplest and well-known quasi-coherent approximation is valid, which reduces by one the effective number of 
dimensions in the transfer equations. The second zone requires the mildly-relativistic effects to be taken into account as they 
determine the frequency redistribution of cyclotron radiation out of the line’s core. In this case a system of two equations for 
the ordinary and extraordinary mode intensities can still be used, and we derive the proper coefficients in these equations. The 
third zone calls for the most comprehensive analysis, with both the mildly-relativistic and QED effects taken into account. In 
this case, the linear mode coupling necessitates to use the complete set of the four transfer equations for the Stokes parameters. 
However, even in a general case, the latter equations have to be used only in the region of linear mode coupling, i.e., close to 
the point of vacuum resonance, and outside this region the equations for incoherent X-, O-modes are valid. 

The complexity of the radiation transfer modeling in the third zone is due to an inhomogeneity of plasma in a region of 
the ’vacuum resonance’, where the modes change their appearance from the plasma-determined polarization to the vacuum- 
determined polarization. For the neutron star atmospheres and at the frequencies near the cyclotron resonance, the transition 
region can be located at the optical depths larger than unity, even for the relatively weak magnetic fields on the order of 
10 11 G. The effect of the linear mode coupling can lead to the observable features in the emergent spectra, though its detailed 
analysis is a subject of future work and will be presented in a follow-up paper. 

In the present paper, we derive the coefficients in the transfer equations both for the X-, O-mode intensities and for the 
Stokes parameters. The coefficients in this form are necessary and sufficient for solving the radiation transfer problem for 
the corresponding regions in the plasma parameter space. The derived equations take into account the effects of the vacuum 
polarization, inclined magnetic field, and relativistic corrections to the cyclotron resonance condition. A detailed analysis of 
these effects, both in the isolated neutron stars and accreting X-ray pulsars, is hindered by a high complexity of the problem. 
However, a progress in observations and a rapid growth of the available computing resources will allow one to use the results, 
obtained in this paper, for building realistic models of the atmospheres of such objects in the near future. 

We also note that all physical effects analysed in this paper can be important for the atmospheres of magnetized white 
dwarfs too, though the influence of the vacuum polarization and the subsequent linear mode coupling is significant only for 
the extreme white dwarfs with a magnetic fields exceeding 10 9 G, which are not observed yet. 

For those readers, who are interested in the numerical analysis of the dispersion properties of the mildly-relativistic 
magnetized plasma, we publish the MATLAB code, which evaluates the components of the permittivity tensor and refraction 
indices. The code can be found at http://www.mathworks.com/matlabcentral/fileexchange/52024 . 
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